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Ultrahigh-energy neutrino-nucleon deep-inelastic scattering and the Froissart bound 
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We present a simple formula for the total cross section cr"^ of neutral- and charged-current deep- 
inelastic scattering of ultrahigh-energy neutrinos on isoscalar nuclear targets, which is proportional 



to the structure function F2 {My/s, 
square of the center-of-mass energy. 



My), where Mv is the intermediate-boson mass and s is the 
The coefficient in the front of {x,Q'^) depends on the 
asymptotic low- a; behavior of F^^. It contains an additional Ins term if scales with a power 
of ln(l/a;). Hence, an asymptotic low-a; behavior ^2'^ oc ln'^(l/a;), which is frequently assumed in 
the literature, already leads to a violation of the Froissart bound on a"^ . 

PACS numbers: 13.15.-l-g, 13.85.Hd, 25.30.Pt, 95.85.Ry 



For more than a decade, large experiments have been 
searching for ultrahigh-energy (UHE) cosmic neutrinos 
[v), with energies Ei, > 10^ GeV, using detectors scan- 
ning for events in large volumes of water, ice, the Earth's 
atmosphere, and the lunar regolith [ij. While no clear 
indication of such an event has yet been reported, exper- 
imental bounds on UHE-neutrino fluxes could be estab- 
lished, which, put together, now cover energies way up to 
10^^ GeV and start to constrain scenarios of astrophysi- 
cal interest. Since these limits directly depend on the to- 
tal cross section a^^ {E^) of UHE-neutrino deep-inelastic 
scattering (DIS) off nucleons {N), it is an urgent task to 
provide reliable theoretical predictions for the latter in 
the asymptotic high-i?^ regime, which lies far beyond 
the one explored by laboratory-based vN DIS experi- 
ments and corresponds to asymptotically low values of 
Bjorken's scaling variable x. This requires extrapolation 
over several orders of magnitude in E^,, for which vari- 
ous approaches exist • These are based on successful 
descriptions of the terrestrial data within the framework 
of perturbative QCD and frequently impose the Froissart 
bound [9^] on cr'^^. According to the latter, unitarity and 
analyticity limit the total cross section of a scattering 
process not to grow faster with energy than In^ s. 

In this Letter, we derive a general formula for ct"^ 
that is remarkably concise and correctly accounts for 
the asymptotic high-energy behavior making it perfectly 
suitable for UHE-neutrino phenomenology. It is propor- 
tional to the DIS structure function {x,Q^), which 
has a well-known representation in terms of parton distri- 
bution functions (PDFs) within the parton model (PM) 
of QCD, with X and the typical energy scale Q appropri- 
ately defined in terms of E^ and My {V — W,Z). To be 
on the conservative side, we assume for the time being, 
as in Refs. that the available experimental data on 

DIS allow for an extrapolation to very high (low) values 
of E^ (x) using an appropriate parameterization of i^^^. 
If the latter rises too steeply as x ^ 0, then possible 
new QCD phenomena, such as gluon saturation or re- 
combination, color glass condensates, multiple pomeron 



exchanges, etc., are expected to enter the stage as a cure 
at x values below those currently probed by DIS experi- 
ments (for a review, see Ref. We shall return to this 
issue below, considering two popular models of screening 

[nlEl. 

Specifically, we consider the charged-current (CC) and 
neutral-current (NC) DIS processes. 



v{k)+N{P) 
v{k)+N{P) 



Kfc') 



X, 



(1) 



respectively, where N — {p + n)/2 denotes an isoscalar 
nucleon target of mass M, X collects the unobserved part 
of the final state, and the four-momentum assignments 
are indicated in parentheses, and introduce the familiar 
kinematic variables 



2q-P' 



y 



k-P' 



(2) 



where q = k — k'. In the target rest frame, we have 
s = M(2K + M) and xy = Q'^/{2ME^,). In the kine- 
matic regime of interest here, the inclusive spin-averaged 
double-differential cross sections of processes ^ are, to 
very good approximation, given by [3| 



dx dy 



GlME, ^ 

A, 

27r 



Ml 



r + M2 



K{y)F!^'', 



(3) 



where i = CC,NC, Gp is Fermi's constant, and K{y) = 
2 — 2y + y^. In the so-called wee parton picture appro- 
priate for the low-a; regime 0, we have Kcc ~ 1 ^nd 
K]<!c = 1/2 — a;m + (10/9)a;^, where Xyj — sin^ 6^, with 6^ 
being the weak mixing angle. Using x^, = 0.231 [l^, we 
have Kj^c = 0.328. The contributions due to the struc- 
ture functions and Fg"^ to the r.h.s. of Eq. dJ) are 
negligibly small in our applications; Fj^^ tends to zero 
as rises [l3|, and F^^ essentially refers to valence 
partons, which hardly contribute in the low- a: regime. 

Detailed inspection of the available £N DIS data (see, 
e.g.. Fig. [T]for ep data from HERA I [15]) suggests that. 



2 



in the limit a; — >■ 0, exhibits a singular behavior of 
the form F^^ {x, Q^) ~ x-^F^^{x, Q^), where (5 is a smah 
positive number and diverges less strongly than any 
power of X, i.e., F2^{x,Q'^)/x^^ — > as a; for 
any positive number A. Assuming a symmetric quark 
sea, as is appropriate for the low- 2: regime, we have 
F^^{x,Q^) = {18 / 5) Fi^{x,Q^), so that the low-x be- 
havior of F2^ carries over to F2^ . 

Imposing the lower cut-off Qq on Q^, the total cross 
sections of processes (HJ are evaluated as 



1 



2ME„ 



1 dx d^ar'^ 



j. X dxdy 



, (4) 



where x = Q^/{2ME^). Inserting Eq. (O, Eq. dH) be- 
comes 



i 


c, X 10^ 


5i X 10^ 


aii X 10^ 


a2i X 10^ 


bi X 10^ 





189.4 


10.90 


-8.471 


12.92 


2.689 


1 


1.811 


6.249 


4.190 


0.2473 


11.63 


2 


-0.6054 


-0.3722 


-0.3976 


1.642 


-0.7307 



TABLE I: The values of the fit parameters appearing in 
Eqs. (Ull), (Uni), and (O. 



integral on the r.h.s. of Eq. ([5]) may be put to zero, so 
that 

M(.,QM + .) = M(l + .) = ^i±^^ (11) 

becomes independent of x and Q^. (2) On the other 
hand, ii S 1, then we have 



G 



■2 f2ME, 



dQ^ 



Ml 



r + Ml 



The inner integral on the r.h.s. of Eq. ([S]) can be rewritten 
as the Mellin convolution K{x)®F2'^ {x, Q^). Exploiting 



the low- a; asymptotic form F2 (2;, Q 



'^F^^{x,Q^) 



explained above, this Mellin transform may be repre- 
sented, at small values of in the factorized form 
M(f,g2,l + (5)^2^^ (i,Q2) up to terms of 0{x) 
Here, 



M{x,Q^,l + 6) = 2 
1 



1 



1 



M(l-f5),(6) 



S{x,Q^) F-^{x,Q^) 



5{x,Q^) S 
''^'x,Q^)Jx V 

2 



M{n)^ I dxx"-^K{x)^— - + ^— 

Jq n - 1 n n + I 

Hence, Eq. ^ becomes 



(8) 



„vN (TP \ ^ 



dQ^ 



Ml 



Q^ + Ml 



xM{x,Q^,l + S)F^^{x,Q^). (9) 
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Because the dependence of F2^ {x,Q^) and hence 
M{x, Q^, 1 + S) is only logarithmic, the factor ]Ml/{Q'^ 
Ml)Y essentially fixes the scale = Ml m, so that 
Eq. ^ simplifies to 



ar^(i?,) ^ ^K,MlM{i,Ml, 1 + 5)F^''{i,Ml), 



(10) 

where i = Ml/{2ME^). This is our master formula. 
Further simplification depends on the actual size of 5, 
and we distinguish two cases. (1) If 5 is not too small, 
so that x^ <^ const, then the lower limit x of the inner 



Mix, 1 + 6)^ M{x, Q\ 1) = J—— - X. (12) 

5[x,Q^) 2 

We note that ^ is determined by the asymptotic low- a; 
behavior of F!^^ . For instance, if F2''^(a;, Q^) oc lnP(l/a;) 
for a: -> 0, then l/,5(a;,g2) = ln(l/a;)/(p + 1) [l^- 

Now we apply Eq. (|10p to the three most popular types 
of i^l^ parameterization, namely the standard PM rep- 
resentation implemented with up-to-date proton PDFs 
[isl [isl [igj , a modification of the impressively simplistic 
log-log form proposed by Haidt (H) [20], and the more 
sophisticated form recently introduced by Berger, Block, 
and Tan (BBT) [6]. While the dependence of the 
PM representation of F2^ is governed by the DGLAP 
evolution, those of the heuristic H and BBT forms are 
directly determined by global fits to experimental data 
covering a wide range. In the low- a: regime, the PM 
parameterization of may be well approximated by 
the following ansatz: 



Fl 



■IN 
2,PM 



(^,Q') = CpM(Q')a;-*™(«'\ 



with 



Cpm(Q^) 
<5pm(0') 



Co 
5o 



ci lnQ2 
5i lnQ2 



C2 In' Q2 
^2 In^ Q2 



(13) 



(14) 
(15) 



where it is understood that is taken in units of 
GeV'. To suppress higher- twist effects, impose the cut 
Q2 > 3.5 GeV^ Fitting Eqs. (IT3l)-([15l) to the result for 
^2^ evaluated at next-to- leading order (NLO) with the 
HERAPDFl.O [13 proton PDFs, we obtain the values of 
Ci and 8i collected in Table ID From Eq. ([T5|) and Table U 
we glean that 



(5pm(A^ 



{Ml,) « 0.37, 



(16) 



so that Eq. is to be used with Eq. jTT]). Us- 

ing the MSTW PI and CTIO ^ PDFs, we obtain 
5^yi{Ml) sa 0.35 and 0.38, respectively. The result- 



ing high-£;i, behavior (J^^ (E^\ol j-<5pm(m^, 
agreement with other studies 



is in good 



For the reader's convenience, we recollect here the BBT 
parameterization of i^2^ appropriate for the range x < 
0.11 [1] relevant for our applications [2l|. It reads 



IN 

2. BBT 



(a;,Q2) ^ (1^2;) 



Ao + Ai(g2)ln 



XP{1 



x{\ ~ Xp) 



A2(Q')ln 



2 Xp{\ - x) 



x{l — Xp) 



(17) 



where = Fp/{1 - xp), with Fp = 0.413 3, and 
A,(g2) = a,o + aa In Q2 + a,2 In' (* = 1,2), (18) 

with the values of listed in Table |T] [8]. Here, Eq. ([TU]) 
is to be used with Eq. ([T2|). and wc find 



1 



JBBT 



{x,Q^) 



ELo Ai\n{xp/x) 



1 , 

- In — . 

3 X 

(19) 



From Eqs. (IT7|) and (IT9t . we glean that, in the high- 



energy limit s 



^2, BBT 



{x,My) oc In s and 



1/Sbbt{S:, Aly) oc Ins. This leads us to the important 
observation that crgB't,^'- ^^ which manifestly violates 
the Froissart bound in contrast to what is stated in 
Refs. This violation of the Froissart bound is at- 

tributed to the presence of the In' x term in Eq. ([TT]) . 
On the other hand, if just rises linearly in Inx as 
X 0, then tr^^ oc In' s in accordance with the Froissart 
bound. 

In fact, this is the case for the original H ansatz (20| . 
B ln(a;o/a;) ln(l + Q^/Qq), which contains just three fit 
parameters. To enable the fit quality to be improved, we 
introduce three more parameters by writing 

F2^^(x,g') = i3o + i?i(Q')ln^, 



i3i(Q') = E&an^(l + ^) 



(20) 



Eq. ([TU]) is again to be used with Eq. p^ . and we obtain 
1 ^XloBJ^'^\xo/x)/{t + l) 



SYi{x,i 



iln^, 
2 X 



(21) 

so that CT^^ oc In' s as it should. Fitting Eq. (gOl) to the 
recent combination of the complete HI and ZEUS data 
sets on i^l^ from HERA I [l^j with the cuts x < 0.01 
and g' > 1.5 GeV', we obtain xq = 0.05791, Ql = 
2.578 GeV', Bq = 0.1697, and the values of 6, hsted in 
TableU with x'/d.o.f. = 422/175 « 2.41. 

Looking at Fig. [U we observe that our fit (solid lines) 
indeed yields a surprisingly good description of the ex- 
perimental data over the full x and g' ranges consid- 
ered. The approximation works particularly well for low 
X and large values and is thus likely to allow for a 
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FIG. 1: Measurement of [T^ compared with our fit using 
the improved H ansatz as well as the PM and BBT results. 



reliable extrapolation to the x and g' ranges relevant 
for UHE-neutrino physics. In fact, switching to the cut 
X < 10^^ reduces x'/d.o.f. by roughly a factor of three, 
to x'/d.o.f. = 58/69 ~ 0.84. For comparison, also the 
PM results evaluated from Eqs. (fT3|) -([T5 l) (dashed lines) 
and the BBT results evaluated from Eqs. (H?]) and (0 
(dash-dotted line) , which are hardly distinguishable from 
one another, are shown in Fig. [TJ 

We now consider i/N DIS with UHE neutrinos. For 
the sake of brevity, we focus our attention on CC DIS. 
The corresponding NC results may be obtained by sub- 
stituting Kcc K]<sc and Mw Mz in our formulas. 
In order to determine the range of validity of our master 
formula (fTO)) for cr^^, we compare it with the exact for- 
mula ([5|) , which requires two-dimensional numerical inte- 
gration, for the PM, BBT, and H cases considered above. 
In each case, we find excellent agreement for E^, values of 
order 10^ GeV and above, which corresponds to x values 
of order 10"'^ and below in This is illustrated for 

the BBT case in Fig. [21 where the approximate evaluation 
of Eq. ([ini) with Eqs. ([T2|) and ([TU) is compared with the 
exact one of Eq. ^ with Eqs. ([T7|) and ([T8|) (dashed line). 
The large- approximation may be somewhat improved 
by evaluating S{x,Q^) hy one-dimensional integration via 
Eq. ^ instead of using Eq. (dotted line). 

The PM and H results for ctqq evaluated from our mas- 



ter formula ([TOj), with Eqs. ([TT]) and ([T6|) in the PM case 
and with Eqs. (IT^ and ([?T|) in the H case, are also dis- 
played in Fig. [21 Comparing them with the corresponding 
BBT result, we observe that all three predictions agree 
relatively well in the range 10'' GeV ^ < 10^ GeV, 
where the high-£^^ approximation is already working 
and the respective i^2^ parameterizations are still con- 



4 




[GeV] 

FIG. 2: Predictions for cfcc{Ev) evaluated from the PM, 
BBT, and H parameterizations of F|^(x,(5^). In the BBT 
case, also the improved high-_E^ approximation and the exact 
evaluation are shown. 



strained by the HERA data. However, these three predic- 
tions steadily diverge as further increases until they 
differ by 1-2 orders of magnitude at typical UHE val- 
ues of El,, reflecting the different low- a; behaviors of the 
respective parameterizations of 

In summary, we derived a novel concise relationship, 
given by Eqs. ([T0|) - ((T2]) . between the total cross section 
of CC and NC vN DIS in the high-£;^ limit and 
the structure function (x,Q'^) in the low- a; limit. It 
is particularly useful for applications to UHE-neutrino 
physics providing reliable predictions in a very quick and 
convenient way as it does. Being given in terms of a 
closed analytic formula, it also allows one to unambigu- 
ously determine if erf ^ resulting from a given functional 
form of satisfies the Froissart bound Q or not, while 
this is hardly possible using the numerical solution of the 
two-dimensional integral in Eq. ([SJ. Specifically, if 
exhibits a low- a; behavior cx ln^(l/x), which corresponds 
to a high-s behavior oc In^s in Eq. (|10l) . then the coef- 
ficient M in that equation produces an additional factor 
oc Ins, so that the Froissart bound is violated for p > 1. 
In fact, this is the case for the BBT 0-11] parameteriza- 
tion of i^2^, for which p ^ 2. On the other hand, the 



H [1^ one is characterized hy p = 1, so that the Frois- 
sart bound is satisfied. This motivated us to update the 
analysis of Ref. [l^l by fitting our improved ansatz ([20|) 
to the recent combination of the complete HI and ZEUS 
data on from HERA I [15]. 

For completeness, we also performed a fit to the PM 
result for evaluated at NLO with an up-to-date set 
of proton PDFs, namely the HERAPDFl.O one obtained 
by HI and ZEUS by fitting their own data |15j|, and pre- 
sented the resulting prediction for a'^^ . As expected, 
the low- a; behavior of the PM result for Fj''^ is too sin- 
gular for crf^ to satisfy the Froissart bound. It is likely 
that the inclusion of nonlinear terms, such as screening 
corrections generated by gluon saturation or recombina- 
tion, in the evolution equations will cure this problem 
10|. In fact, considering the Ayala-Gay-Ducati-Levin 
111 and the generalized Golec-Biernat-Wiisthoff (GBW) 
1^ models of saturation, where, due to their specific 
gluon densities, cx ln^(l/a:) withp = 1 andp — 0, 
respectively, we obtain a^^ oc In^"*"^ s, where the second 
additional logarithm arises from the integration in 
Eq. Thus, saturation strongly modifies the power- 
like perturbative asymptotics of total cross sections and 
has the potential to restore the Froissart bound, as in the 
case of the GWB model. 

Future measurements of lyN DIS with UHE neutrinos 
will eventually provide direct access to the low- a; asymp- 
totic behavior of far beyond the reach of accelera- 
tor experiments, and our new relationship will provide a 
convenient tool to expose it. On the theoretical side, one 
important lesson to be learned from our specific exam- 
ple, where total cross sections could be simply related to 
structure functions in the framework of perturbation the- 
ory, is that the direct application of the Froissart bound 
to structure functions represents a potential pitfall, of 
which we wish to caution the reader. 
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